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Equilibrium Distribution of Heavy Quarks in Fokker- Planck Dynamics 
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We obtain within Fokker-Planck dynamics an explicit generalization of Einstein's relation between 
drag, diffusion and equilibrium distribution for a spatially homogeneous system, considering both 
the transverse and longitudinal diffusion for dimension n > 1. We then provide a complete char- 
acterization of when the equilibrium distribution becomes a Boltzmann/Jiittner distribution, and 
when it satisfies the more general Tsallis distribution. We apply this analysis to recent calculations 
of drag and diffusion of a charm quark in a thermal plasma, and show that only a Tsallis distri- 
bution describes the equilibrium distribution well. We also provide a practical recipe applicable to 
highly relativistic plasmas, for determining both diffusion coeflicients so that a specific equilibrium 
distribution will arise for a given drag coefficient. 

PACS number(s): 12.38.Mh, 05.10.Gg, 52.65.Ff, 05.20.Dd 



The velocity distribution of objects subject to a (ther- 
mal) background plays an important role in a number 
of scientific fields, including plasma physics and astro- 
physics 10], nuclear physics |^,^, and, more generally, in 
kinetic theory [0-^. The Fokker-Planck equation is a 
popular tool to study this distribution. It can be mo- 
tivated in a number of ways. One method is to create 
a Langevin equation l^,^, which describes the stochas- 
tic behavior of a single object propagating with ran- 
dom noise. Another method comes by taking a mas- 
ter equation, such as the linearized Boltzmann-Vlasov 
equation and performing a Landau soft-scattering 
approximation. Finding recent theoretical calculations 
for the transport coefficients of the Fokker-Planck equa- 
tion based on a microscopic theory [^jsj, we recognized 
the need to establish a simple procedure for understand- 
ing the relation between the transport coefficients of the 
Fokker-Planck equations as determined in such micro- 
scopic calculations, and the resulting properties of the 
equilibrium distribution. 

After a brief summary of the recent developments in 
the Fokker-Planck studies of equilibrating heavy quarks 
in a quark-gluon plasma, we generalize and apply a well- 
known relation between kinetic coefficients and the equi- 
librium distribution of the Fokker-Planck equation in a 
spatially homogeneous environment In that way we 
are able to relate the drag and diffusion coefficients to 
the shape of the equilibrium distribution. We stress the 
importance of including both transverse and longitudinal 
diffusion to maintain a consistent equation. A simple test 
follows which exactly determines when the equilibrium 
distribution obeys Boltzmann/Jiittner statistics or the 
more general Tsalfis statistics We also discuss how 
to choose the transport coefficients in order to attain the 
Boltzmann/Jiittner distribution, and address some issues 
related to the difference between the stopping power and 
the drag and diffusion coefficients. 

The statistical properties of an ensemble of objects 
(particles) can be expressed in terms of the one-particle 



distribution function, f{x,p,t). This density, when mul- 
tiplied by the 2n-dimensional phase-space volume ele- 
ment (f^x d"p, gives the probability of finding the object 
in this infinitesimal region of phase-space. We have intro- 
duced the dimensionality n explicitly, and we will primar- 
ily pursue the physical case n = 3, with the case n = 1 
also of interest due to its exceptional character. We as- 
sume that f{x,p,t) obeys a Boltzmann-Vlasov master 
equation of the form |^,^: 



d"k[W{p + k,k)f{x,p + k,t) 
-W{p,k)f{x,p,t)]. (1) 



where: 

/ = f{x,P,t), 



dx p u dp ^ 

In the non-relativistic limit, E m, but otherwise our 
notation is applicable to both classical and relativistic 
mechanics. The collision term has two parts: in the first 
gain term the transition rate W{pi , k) represents the rate 
that a particle with momentum pi = p + k loses momen- 
tum k due to reactions with the background. The second 
term represents loss due to scattering out. The collision 
term is strictly local, depending only on momenta of par- 
ticles, but it depends on position x indirectly, because W 
incorporates any background inhomogeneity. 

Expanding the gain term about p to second order in k 
leads to the Fokker-Planck equation: 

a/ 

dt 



df df d ^ ^ d d ^ ^ 
dxt dpi dpi dpi dpj 



where we are using the Einstein summation convention 
for repeated indices i and j. We have introduced the 
transport coefficients of drag and diffusion, respectively: 



(4) 



A,^A,{p) = J d''khW{p,k), 
= B,^{p) ^Ifd^k hkj W{p, k) 
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It is generally believed that the Fokker-Planck equa- 
tion describes well the approach to (thermal) equilib- 
rium. However, since we shall find that this is not guaran- 
teed, we record yet another independent way to motivate 
the form of the Fokker-Planck Eq. (|^), the Ito-Langevin 
method. Consider the Langevin system of equations: 



dx 



P_ 

E 



dt 



= F, (.f) + G, {x, p) + {x, p)Vj (t) , (5) 



where the noise term ff is Gaussian white noise with 
{rii{t)) = 0, and {Viit)vjit')) = S{t - t'). Using Ito's 
formula one shows that this Langevin system corresponds 
to the Fokker-Planck equation M : 



dl 
dt 



Pi df 
E dxi 



F, 



df -dGif 



dpi 



dpi 



d\lDD^),,f 
dpidpj 



(6) 



—Gi and Bij <-> 



Thus, we see that we can identify A, 

While the computation of Dij is not obvious in the 
Langevin formulation, the master equation approach 
gives precise formulas. Written in terms of the two body 
collision reaction matrix elements the drag and dif- 
fusion, according to Eq. (0) , are : 



1 f d^k 
"^'^P^ - 2eJ J2^. 



d^k' 



d^p' 



f2EkJ {2nf2Ek'J {2nf 2Ep, 
- J2 \M\\2n)^S\p + k-p'- k')[p, p'MMk) (7) 



^{{Pi-p[)); 
Bijip) = i^{{{Pi - p'i){P] -P'j))) 



(8) 
(9) 



In our case, the incoming particle is different from the 
background. For each background species, there is a 
similar additive contribution to the collision integral in 
Eq. (0). Moreover, as long as the background does not 
distinguish discrete quantum numbers of incoming par- 
ticles (such as spin), we can combine these properties 
in one distribution /. In such a case, we need to aver- 
age over the initial reaction channels and sum over the 
open final channels, akin to the evaluation of the cross 
section, hence the degeneracy factor 7"^ of the fore- 
ground particle. In Eq. (|^) g{k) is the particle density 
of the background, assuming that there is a single type 
of particle, and g{k) = [1 ± g{k)] represents a Bose en- 
hancement/Pauli suppression factor for scattered back- 
ground particles, as appropriate. We assume that the 
background has equilibrated at the temperature Tb. 

We are now prepared to study the equilibrium distri- 
bution and its relation with the drag and diffusion coeffi- 
cients. We study the simplest possible case of a spatially 
homogeneous distribution. In absence of vectors other 
than p the values of Ai and Bij , which depend function- 
ally on p and the background temperature T, must be 
of the form, where p^ = pf (summation convention is 
always implied) : 



A,{p,T) = PiA{p,T) , 
B,,{p,T)=(5,,-^]B^{p,T) 



(10) 



'^B,{p,T). (11) 



is the longitudinal and B± the transverse diffusion 
coefficient. In terms of microscopic reaction amplitudes, 
these three functions are defined by the following expres- 



Aip) = ((!)) - 

P 



B±ip) 
Biiip) 



m-p'fi) 



2{{p-p'))+p'{{l)) 



(12) 
(13) 
(14) 



With no external forces and a homogeneous back- 
ground, Eq. (jsl) reads: 



(15) 



A natural requirement for /oq, detailed balance, is for the 
probability current P to vanish |^: 

A,{p,T) ^ Bij{p,T)— i , (16 

dPo dpj 

where we wrote the equilibrium distribution as: 



/oq(p;r,g) =iVexp(-$(p;r,g)) 



(17) 



where T, q are parameters that may be needed to char- 
acterize the shape of the distribution. R. Graham and 
H. Haken ||l^,|ll| provide a more general approach which 
classifies when such a simplification is valid, and how to 
extend the condition Eq. ([l^) when it is not valid. For 
our case, we have verified that Eq. (16) is valid. 

Using Eqs. (|l^,|ll|) and the fact that in the spatially 
homogeneous case the equilibrium distribution depends 
only on p = \p\, Eq. ( p^ ) becomes 

. ld«> , , IdBn 
A{p,T)^-—B\\{p,T)--—^ 
p dp p dp 

n — 1 

-^iB\\iP,T)-B^ip,T)), (18) 

which is the desired relation between the shape of equi- 
librium distribution and the three drag/diffusion coeffi- 
cients. The special case considered by Einstein arises in 
the classical problem of a particle that travels through an 
ideal heat bath and undergoes linear damping (Rayleigh's 
particle). Substituting the coefhcients A = and B± — 
= D into the relation (p^, we obtain the Boltzmann 
equilibrium distribution Eq. ( p^ with $ = p^ /2mkT 
only if Einstein's well-known drag-diffusion relation for 
Brownian motion 7 = D/mkT is satisfied. 
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When the equihbrium distribution is known a priori 
and if the diffusion coefficients are also known, then it is 
an easy matter to use Eq. (|f^) to find the unique, con- 
sistent drag coefficient. However, the reverse is not true 
for n > 1: given the equihbrium distribution and the 
drag coefficient, there are two diffusion coefficients which 
must be simultaneously determined, which is in general 
not possible without a further assumption. One 'popu- 
lar' method to do this is to assume that the tensor Bij 
is diagonal, that is B±^ = B\\ and then to solve the linear 
first order differential equation to obtain Bn : 



(19) 



An alternative method, motivated when Tb is the only 
dominant energy scale, is to calculate the longitudinal 
diffusion term B\\ in terms of the drag and then to deter- 
mine the transverse diffusion term B±^ : 



Bn 



pA 

d^/dp 



Bi 



Bn 



P 



dB 



n — 1 dp 



(20) 



We note that when = const, both recipes give the 
same result. These approximate relationships are suffi- 
cient to guarantee correct equilibrium behavior, but may 
not represent accurately the dynamical evolution, and 
one should always require non-negative diffusion. 

While discussing the relationships between the drag 
and diffusion coefficients of the Fokker-Planck equation, 
we note another interesting relation between these coeffi- 
cients and the frequently discussed stopping power. The 
stopping power measures the energy loss per unit dis- 
tance traveled, and is equal to the energy loss per unit 
time divided by the particle speed. Thus, in terms of the 
elementary matrix elements [Q: 

-^ = ((^)) = ((^— ^)). (21) 
dx V p 

When combined appropriately with Ai , a relativistically 
invariant (scalar) quantity is found: 

dF, 

p.A+p^ = {{p- ip-p') - E' + EE')) (22) 
dx 

= -l{{{p,-p'n')), 



where we use four- vector notation. Eq. (22) shows that 
the stopping power and the drag coefficient are, in gen- 
eral, two independent quantities. To connect them we 
need to evaluate also: 



Boo ^ {{{E - E'f)) , ^ {{E - E')) . 



(23) 



In the non-relativistic limit these two new quantities are 
relatively small. The energy loss can be expressed as: 



dE _ piAj 
dx 



B: 



00 



Bi; 



dE 



(24) 



For the Rayleigh particle considered earlier, noting 
Bqo — > 0, this corresponds to an energy loss: 



dE 



dx p 



2m I p" 
2m 



3kT\ 



7m(f 



(25) 



which vanishes precisely for a thermal velocity. 

The problem that we are facing even in the simple spa- 
tially homogeneous case is that the Fokker-Planck coeffi- 
cients cannot simply be chosen to assure that the 'correct' 
equilibrium distribution results but have already been 
obtained in terms of elementary collision reaction am- 
plitudes, see Eqs. (p^-p^), and thus the resulting equi- 
librium distribution is fixed, as can be seen solving and 
integrating Eq. ( p^ to obtain <I>. Since the drag and dif- 
fusion coefficients are not evaluated exactly but in some 
valid approximation, typically applying a perturbative 
expansion, it is more appropriate to analyze the result- 
ing distribution in terms of some useful class. We con- 
sider the class of Tsallis statistics [|j , which depends on a 
temperature-like quantity T and on a parameter q, which 
measures the degree of extensivity of entropy in the sys- 
tem: 



/eq-7V[i-(i-g)i?(p)/r]- 

d<P dE 1 

~ 'd^T- {l-q)E ' 



(26) 
(27) 



The Boltzmann distribution arises when q ^ 1. Substi- 
tuting into Eq. dlS), we obtain: 



T+iq-l)E 



dE 



B 



dp pA- 



dp 



(28) 



Whenever the ratio given by the right-hand side of 
Eq. ( p8| ) becomes linear in E, then Tsallis statistics de- 
scribe the stationary distribution. When the ratio is 
constant, then a Boltzmann/ Jiittner distribution suf- 
fices. We note that for the special case n — 1 and non- 
relativistic dynamics dE /dp — v, Eq. ( p8[ ) was obtained 
recently within the Langevin formulation of the Fokker- 
Planck dynamics jl^ . 

We consider now the drag and diffusion coefficients for 
a charm quark with mass TOc = 1-5 GeV interacting with 
thermal gluons at Tb = 500 MeV calculated using per- 
turbative QCD techniques [||J|]. We have gone to con- 
siderable length to assure that these results apply 0. 
Diamonds in Fig.[| show the ratio Eq. (p8|). The linear 
regression fit (straight line) shows that the parameters 
best describing the distribution as a Tsallis distribution 
are q = 1.114 and Tt = 135.2 MeV. The dashed horizon- 
tal line in Fig. |l] corresponds to the Boltzmann/Jiittner 
distribution {q — 1) and Tt = Tb), which we were ex- 
pecting to find. The difference to the actual distribution 
appears to be significant. 
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3.5 
E [GeV] 

FIG. 1. Calculated data (diamonds) and linear fit for the 
ratio in Eq. ( psj ) for a charmed quark ttIc = 1.5 GeV thermal- 
izing in gluon background at Tb — 500 MeV. Dashed line: re- 
sult expected for a Boltzmann/Jiittner distribution, T = Tb. 

The more practical question is what the charmed quark 
spectrum would actually look like. This is shown in Fig. || 
where a solid line shows the Tsallis distribution as ob- 
tained above, compared to Boltzmann/Jiittner shape for 
Tb = 500 MeV and = 1.5 GeV. Assuming that the 
Tsallis shape would be measured, the spectrum would 
reveal two components: at low E a 'cold' Boltzmann dis- 
tribution, and for high E a power-law with /cq oc E", 
where in our case n — — —8.8. 




3.0 
p [GeV] 

FIG. 2. Normalized equilibrium spectra: Tsallis distribu- 
tion (solid line) and Boltzmann/Jiittner distribution (dashed 
line) at Tb = 500 MeV, with = 1.5 GeV. 

In summary, we developed tools which allow one to 
identify within Fokker-Planck dynamics the equilibrium 
distribution for given (calculated) drag and diffusion co- 



efficients, or when the stationary distribution is known, 
to determine the drag or, as a recipe for n > 1 di- 
mensions, both longitudinal and transverse diffusion co- 
efficients. We have then shown that thermalization of 
charmed quarks in a quark-gluon plasma leads to the 
two parameter Tsallis distribution, and have determined 
the pertinent parameters for the published microscopic 
drag/diffusion coefficients. It is at present unclear if the 
resilience of the Tsallis statistics is a fundamental fea- 
ture of the Fokker-Planck dynamics in background rela- 
tivistic plasma p5| . It is important to note that only a 
major change in the transport coefficients from the re- 
sults of the microscopic calculations will lead to a Boltz- 
mann/Jiittner equilibrium distribution. 
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